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We introduce the basic techniques used for the analysis of three dimensional and 
two dimensional galaxy samples. We report the correlation analysis of various red- 
shift surveys which shows that the available data are consistent with each other 
and manifest fractal correlations (with dimension _D ~ 2) up to the present obser- 
vational limits without any tendency towards homogenization. This result points 
to a new interpretation of the galaxy number counts. We show that an analysis of 
the small scale fluctuations allows us to reconcile the correlation analysis and the 
number counts in a new perspective which has a number of important implications. 



1 Introduction 

In this lecture we briefly introduce the basics concepts of fractal geometry and 
the methods of correlation analysis, that are usually used in Statistical Me- 
chanics. First we introduce the methods to compute the full correlation func- 
tion in three dimensional samples characterized by a wide solid angle, then 
we consider the determination of the space density in very deep and narrow 
surveys. Such an analysis allows us to study and characterize the effect of the 
small scale fluctuations in the determination of non-average quantities, as the 
radial density. Having clarified these effects we move to the interpretation of 
the galaxy number counts, i.e. the number of galaxies with a certain apparent 
magnitude versus the magnitude itself. These counts must be examined with 
great care, because, also in this case, it is not possible to perform an average 
over different observers. We show that an analysis of the small scale fluctua- 
tions allows us to reconcile the correlation analysis and the number counts in 
a new perspective which has a number of important implications. 

Finally we summarize our main conclusions. For a more detailed discussion 
we refer the reader to Coleman & Pietronero(1992) and Baryshev et aL(1994) 
for a basic introduction to this approach, and to Pietronero et o/.(1997) and 
Sylos Labini et aL(1997) for a review on the more recent results (see Davis 
(1997) for a different point of view on this subject). 
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2 Space distribution 



We first introduce some basic definitions. If L is tlie absolute or intrinsic 
luminosity of a galaxy at distance r, this appears with an apparent flux 

For historical reasons the apparent magnitude m of an object with incoming 
flux / is 

TO = —2.5 logj^Q / + constant, (2) 

while the absolute magnitude M is instead related to its intrinsic luminosity 
Lby 

M — —2.5 logj^o ^ + constant'. (3) 

From Eq.0 it follows that the difference between the apparent and the abso- 
lute magnitudes of an object at distance r is (at relatively small distances, 
neglecting relativistic effects) 

TO - M = 5 logio r + 25 (4) 

where r is expressed in Megaparsec {IMpc — 3 • 10^^ cm). 

A catalog is usually obtained by measuring the redshifts of the all galaxies 
with apparent magnitude brighter than a certain apparent magnitude limit 
TO/im, in a certain region of the sky defined by a solid angle Q. An important 
selection effect exists, in that at every distance in the apparent magnitude lim- 
ited survey, there is a definite limit in intrinsic luminosity which is the absolute 
magnitude of the fainter galaxy which can be seen at that distance. Hence at 
large distances, intrinsically faint objects are not observed whereas at smaller 
distances they are observed. In order to analyze the statistical properties of 
galaxy distribution, a catalog which does not suffer for this selection effect 
must be used. In general, it exists a very well known procedure to obtain a 
sample that is not biased by this luminosity selection effect: this is the so- 
called "volume limited" (VL) sample. A VL sample contains every galaxy in 
the volume which is more luminous than a certain limit, so that in such a sam- 

Ilfi there is no incompleteness for an observational luminosity selection effect 
u. Such a sample is defined by a certain maximum distance Rvl and the 
absolute magnitude limit Mvl given by 

MvL = rnitm - 5 log^Q Rvl - 25 - A{z) (5) 

where A{z) takes into account various corrections (K-corrections, absorption, 
relativistic effects, etc.), and mum is the survey apparent magnitude limit. 



2 



Different VL samples, extracted from one catalog, have different Rvl, and 
deeper is a VL sample, larger its Ry l ■ 

The measured velocities of the galaxies have been expressed in the pre- 
ferred frame of the Cosmic Microwave Background Radiation (CMBR), i.e. 
the heliocentric velocities of the galaxies have been corrected for the solar 
motion with respect to the CMBR, according with the formula 

V — Vm + 316cos6 kms~^ (6) 

where v is the corrected velocity, Vm is the observed velocity and is the 
angle between the observed velocity and the direction of the CMBR dipole 
anisotropy [a — 169.5° and 5 = —7.5°). From these corrected velocities, we 
have calculated the comoving distances r{z), with for example go = 0-5, by 
using the Mattig's relation □ 

r{z) = 6000 I 1 , ^ ) h~^Mpc . (7) 




In general, we have checked that the results of our analysis depend very weakly 
on the particular value of go adopted, except very deep surveys, and we have 
also used the simple linear relation 

r - cz/Ho . (8) 

In the nearby catalogs there is no any sensible change by using Eq.H instead 
of Eq|. Hereafter for the Hubble constant we use the value Hq = 100 • /i ■ 
kmsec~^ Mpc~^ . 

All the analyses presented here have been performed in redshift space and 
we have not adopted any correction to take into account the eventual effect of 
peculiar velocities (local distortion to the Hubble flow). However we point out 
that peculiar velocities have an amplitude up to ^ 500 lOOO/cmsec^^ and 
then their effect can be important only up to 5 lOh^^Adpc, and not more 
(see the lecture of A. Szalcy in these Proceedings) 

We briefly mention the characteristics of the galaxy luminosity distribution 
that is useful in what follows analyses. The basic assumption we use to compute 
all the following relations is that: 

y{L,r) = c^{L)p{r) , (9) 

i.e. that the number of galaxies for unit luminosity and volume j^(L, r) can be 
expressed as the product of the space density p(f) and the luminosity function 
4>{L) (L is the intrinsic luminosity). This is a crude approximation in view 
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of the multifractal properties of the distribution (correlation between position 
and luminosity) Q. However, for the purpose of the present discussion, the 
approximation of Eq.^ is rather good and the explicit consideration of the 
multifractal properties have a minor effect on the properties we discuss Q. 

To each VL sample (limited by the absolute magnitude Mvl) we can 
associate the luminosity factor 

^{Mvl) = / (l){M)dM (10) 



that gives the fraction of galaxies present in the sample. Hereafter we adopt 
the following normalization for the luminosity function 

$(oo) = / (t){M)dM = 1 (11) 

J —oo 

where Mmm ~ —10-; — 12 is the fainter galaxy present in the available samples. 
The luminosity factor of Eq.|l^ is useful to normalize the space density in dif- 
ferent VL samples which have different Mvl (Eq||). The luminDsity function 
measured in real catalogs has the so-called Schecther like shape eI 

0(M)dM = A ■ 10-°-'*(*'+i)^'^e-i"''"'"'""'dM (12) 

where S « —1.1 and M* « — lO.si, 0, and the constant A is given by the 
normalization condition of Eq.|ll|. 

Eq.| and eqjl^ are useful for the normalization of the average (conditional) 
density in different VL samples. In fact, in view of eq.|[ the difference of the 
average density in different VL samples, is only due to the luminosity factor 
given by eq.p^. 

2.1 Full correlation analysis 

In this section we mention the essential properties of fractal structures because 
they will be necessary for the correct interpretation of the statistical analysis. 
However in no way these properties are assumed or used in the analysis itself. 

A fractal consists of a system in which more and more structures appear 
at smaller and smaller scales and the structures at small scales are similar to 
the ones at large scales. The first quantitative description of these forms is the 
metric dimension. One way to determine it, is the computation of mass-length 
relation. Starting from an point occupied by an object of the distribution, 
we count how many objects A^(r\_(j"mass") are present, in average, within a 



volume of linear size r (" length" ) I 

<N{r)>^B-r° (13) 
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D is the fractal dimension and characterises in a quantitative way how the 
system fills the space. The prefactor B depends to the lower cut-offs of the 
distribution; these are related to the smallest scale above which the system is 
self-similar and below which the self similarity is no more satisfied. In general 
we can write: 

(14) 

where is this smallest scale and N^, is the number of object up to r*. For 
a deterministic fractal this relation is exact, while for a stochastic one it is 
satisfied in an average sense. Eq.(|l^) corresponds to a average behaviour of 
N{r), that is a very fluctuating function; a fractal is, in fact, characterised 
by large fluctuations and clustering at all scales. We stress that eq.(|l^) is 
completely general, i.e. it holds also for an homogeneous distribution, for 
which D — 3. From eq.(|l3|), we can compute the average density < 71 > for a 
sample of radius Rs which contains a portion of the structure with dimension 
D. Assuming for simplicity a spherical volume ( V{Rs) — (4/3)7ri?g), we have 

If the distribution is homogeneous {D = 3) the average density is constant 
and independent from the sample volume; in the case of a fractal, the average 
density depends explicitly on the sample size Rs and it is not a meaningful 
quantity. In particular, for a fractal the average density is a decreasing function 
of the sample size and < n >— > for i?^ — s- 00. 

It is important to note that eq.([l^) holds from every point of the system, 
when considered as the origin. This feature is related to the non-analyticity of 
the distribution. In a fractal every observer is equivalent to-any other one, i.e. 
it holds the property of local isotropy around any observerll3. 

The first quantity able to analyze the spatial properties of point distri- 
butions is the average density. Coleman & Pietronero (1992) proposed the 
conditional density defined as: 

^ <n{r + n)n{n) >. 
< n> 

where the index i means that the average is performed over the points 
of the distribution. In other words, we consider spherical volumes of radius 
r around each points of the sample and we measure the average density of 
points inside them. Such spherical volumes have to be fully contained in the 
sample boundaries. In Eq.^ < n > is the average density of the sample; 
this normalisation does not introduce any bias even if the average density is 
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sample-depth dependent, as in the case of fractal distributions, as one can see 
from Eq. |l^. T{r) (Eq. |l^) can be computed by the following expression 

r('-) = HFi ES" sfe "« + ••')*" = 

where M(r) is the number of points which contribute which are at distance <; r 
from any of the boundaries of the sample. T{r) is a smooth function away from 
the lower and upper cutoffs of the distribution (r* and the dimension of the 
sample). From Eq.(p7|), we can see that r(r) is independent from the sample 
size, depending only by the intrinsic quantities of the distribution {B and D). 
If the sample is homogeneous,!? = 3, r(r) is constant. If the sample is fractal. 



then D < 3, 7 > ai 



,r(r) is a power law. For a more complete discussion 



mc i (r) 

we refer the reader toEl,a. If the distribution is fractal up to a certain distance 
Aq, and then it becomes homogeneous, we have that: 

rW = ^r^-' r < Ao 
An 

r(0 - ^^0^' ^ > ^0 (18) 

It is also very useful to use the conditional average density defined as: 



r*(r) = / ATrr'^T(r')dr' (19) 



This function produce an artificial smoothing of r(r) function, but it correctly 
reproduces global properties El. 

Given a certain spherical sample with solid angle Q and depth i?s, it is 
important to define which is the maximum distance up to which it is possible 
to compute the correlation function (r(r) or ^(r)). As discussed in Sylos 
Labini et aL(1997), we limit our analysis to an effective depth Reff that is of 
the order of the radius of the maximum sphere fully contained in the sample 
volume (Fig.|l|). For example for a catalog with the limits in right ascension 
{ai < a < 02) and declination (5i < S < 62) we have that 

- TT7M69J2) ^ ^ 

where 69 — min{a2 ~ ai, 62 — Si). In such a way we do not consider in the 
statistics the points for which a sphere of radius r is not fully included within 
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'R. 

Figure 1: Upper part: A typical cone diagram for a wide angle galaxy catalog (e.g. CfA, 
SSRS, Perseus-Pisces). The depth is R^- The effective depth is R^ff and it corresponds 
to the radius of the maximum sphere fully contained in the sample volume {Reff ^ Rs)- 
Bottom part: A typical cone diagram for a narrow angle galaxy catalog (e.g. LCRS, ESP). 

In this case R^ff « Rd- 



the sample boundaries. Hence we do not make use of any weighting scheme 
with the advantage of not making any assumption in the treatment of the 
boundaries conditions. For this reason we have a smaller number of points and 
we stop our analysis at a smaller depth than that of other authors. 

The reason why r(r) (or ^(r) see in what follows) cannot be computed for 
r > Reff is essentially the following. When one evaluates the correlation func- 
tion (or power spectrum) beyond Reff, then one makes explicit assumptions 
on what lies beyond the sample's boundary. In fact, even in absence of correc- 
tions for selection effects, one is forced to consider incomplete shells calculating 
r(r) for r > Reff, thereby implicitly assuming that what one does not find 
in the part of the shell not included in the sample is equal to what is inside 
(or other similar weighting schemes). In other words, the standard calculation 
introduces a spurious homogcnization which we are trying to remove. 

If one could reproduce via an analysis that uses weighting schemes, the 
correct properties of the distribution under analysis, it would be not necessary 
to produce wide angle survey, and from a single pencil beam dcc^p survey it 
would be possible to study the entire matter distribution up to very deep scales. 
It is evident that this could not be the case. By the way, we have done a test 
on the homogcnization effects of weighting schemes on artificial distributions 
as well as on real catalogs, finding that the flattening of the conditional density 
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is indeed introduced owing to the weighting, and does not correspond to any 
real feature in the galaxy distribution □. 



2.2 Standard analysis 

At this point it is instructive to consider the behaviour of the standard correla- 
tion function f (r). Coleman & Pietronero (1992) clarify some crucial points of 
the such an analysis, and in particular they discuss the meaning of the so-called 
"correlation length" ro found with the standard approachElEj and defined by 
the relation: 

aro) = 1 (21) 

where 

^ <n(rl)n(rl+rO>. _ ^ 
< n >2 

is the two point correlation function used in the standard analysis. If the 
average density is not a well defined intrinsic property of the system, the 
analysis with ^(r) gives spurious results. In particular, if the system has fractal 
correlations, the average density is simply related to the sample size as shown 
by Eg. (|l5|) . In other words, it is meaningless to define the correlation length 
of the distribution by comparing the average correlation < n{f^i)n(fi + r) >; to 
the average density of the sample < n if the latter depends jpn the sample 
volume. The expression of the ^(r) for a fractal distribution isu 

e(r)-((3-7)/3)(r/i?.)-^-l (23) 

where Rg (the effective sample radius) is the radius of the spherical volume 
where one computes the average density from Eq. (|l5|). From Eq. ( p3| ) it 
follows that: i.) the so-called correlation length rg (defined as $(ro) = 1) is a 
linear function of the sample size Rs 

ro = ii3~j)/6)^Rs (24) 

and hence it is a quantity without any correlation meaning, but it is simply 
related to the sample size. 

ii. ) the amplitude of the ^(r) is: 

A(i?,) = ((3-7)W (25) 

iii. ) ^(r) is a power law only for 

((3-7)/3)(r/i?.)-^»l (26) 
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hence for r < ro: for larger distances there is a clear deviation from a power 
law behavior due to the definition of S,{r). This deviation, however, is just 
due to the size of the observational sample and does not correspond to any 
real change of the correlation properties. It is clear that if one estimates the 
exponent of ^(r) at distances r > tq, one systematically obtains a higher value 
of the correlation exponent due to the break of ^(r) in the log-log plot. This 
is actually the case for the analyses performed so far: in fact, usually, ^(r) 
is fitted with a power law in the range O.Srg ^ r < 2ro, where we get an 
higher value of the correlation exponent. In particular, the usual estimation 
of this exponent by the ^(r) function leads to 7 « 1.7, different from « 1 
(corresponding to « 2), that we found by means of the T{r) analysisQ 

2.3 Normalization of the average density in different surveys 

We have already defined the luminosity factor which is associated to each VL 
sample. As long as the space and the luminosity density can be considered 
independent, the normalization of r(r) in different VL samples can be simply 
done by dividing their amplitudes for the corresponding luminosity factors. 
Of course such a normalization is parametric, because it depends on the two 
parameters of the luminosity function S and M* For a reasonable choice of 
these two parameters we find that the amplitude of the conditional and radial 
density matches quite well in different VL samples. In particular: 

• (1) For CfAl , PP, ESP, LEDA and APM the parameters of the luminosity 
function are 5 = —1.1 and AI* = —19.5. All these surveys have been 
selected in the B-band, and hence the luminosity function is the same. 

• (2) For SSRSl the selection criteria have been chosen in the apparent 
diameter d, and hence we use the linear-diameter function, rather than 
the luminosity function for the normalization of the amplitude in different 
VL samples. We use the diameter function: 

N{D)dD = No ■ exp{-5d ■ D)dD (27) 

where DJis, the absolute diameter (in Kpc), 64 = 0.109 and iVo is a 
constant 0. The shape of this linear-diameter function is similar to the 
Schecther one for magnitudes. 

• (3) Galaxies in LCRS have been measured in the r band, so that we 
have used the luminosity function given by Schetman et aZ.(1996), i.e. a 
Schecther like function with parameters 6 = —0.9 and M* = —20.03. 
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Figure 2: The spatial conditional average density r*(r) computed in some VL samples 
of Perseus-Pisces, CfAl, ESP and LEDA and normalized to the corresponding luminosity 

factor. 



• (4) The two IRAS surveys have been measured in the near infrared. In 
this case we have normahzed the amplitude of the density using the IR 
luminosity function 113. 

We find that in all the cases (l)-(4) the density amplitude in the different 
VL samples match quite well (see Figs.||-^). 

The main data of our correlation analysis are collected in Fig. I (left part) 
in which we report the conditional density as a function of scale for the various 
catalogues. The relative position of the various lines is not arbitrary but it is 
fixed by the luminosity function, a part for the cases of IRAS and SSRSl for 
which this is not possible. The properties derived from different catalogues are 
compatible with each other and show a power law decay for the conditional 
density from \h~^Mpc to 150/i^^A/pc without any tendency towards homog- 
enization (flattening). This implies necessarily that the value of tq (derived 
from the ^ (r) approach) will scale with the sample size Rs as shown also from 
the specific data about rg of the various catalogues. The behaviour observed 
corresponds to a fractal structure with dimension Z) ~ 2. (The data beyond 
^ 200h~^Mpc have been obtained by measuring the radial density and are 
therefore weaker from a statistical point of view - see below) The smaller value 
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Figure 3: The spatial conditional average density r*(r) computed in some VL samples of 
APM, LCRS, SSRSl, IRAS l.2Jy and normalized to the luminosity factor. 



of CfAl was due to its limited size. An homogeneous distribution would cor- 
respond to a flattening of the conditional density which is never observed It is 
remarkable to stress that the amplitudes and the slopes of the different surveys 
match quite well. From this figure we conclude that galaxy correlations show 
very well defined fractal properties in the entire range 0.5 1000h~^Mpc with 
dimension D — 2 ± 0.2. Moreover all the surveys are in agreement with each 
other. 

It is interesting to compare the analysis of Fig.^ with the usual one, made 
with the function ^(r), for the same galaxy catalogs. This is reported in FigjH 
and, from this point of view, the various data the various data appear to be 
in strong disagreement with each other. This is due to the fact that the usual 
analysis looks at the data from the prospective of analyticity and large scale 
homogeneity (within each sample). These properties have never been tested 
and they are not present in the real galaxy distribution so the result is rather 
confusing (Fig.|^). Once the same data are analyzed with a broader perspective 
the situation becomes clear (Fig.||) and the data of different catalogs result in 
agreement with each other. It is important to remark that analyses like those 
of Fig.^ have had a profound influence in the field in various ways: first the 
different catalogues appear in conflict with each other. This has generated the 
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Figure 5: Traditional analyses based on the function ^(r) of the same galaxy catalogs of the 
previous figure. The usual analysis is based on the a priori untested assumptions of analytic- 
ity and homogeneity. These properties are not present in the real galaxy distribution and the 
results appear therefore rather confusing. This lead to the impression that galaxy catalogs 
arc not good enough and to a variety of theoretical problems like the galaxy-cluster mis- 
match, luminosity segregation, linear and non-linear evolution, etc.. This situation changes 
completely and becomes quite clear if one adopts the more general conceptual framework 
that is at the basis the previous figure 
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concept of not fair samples and a strong mutual criticism about the validity 
of the data between different authors. In the other cases the discrepancy 
observed in Fig.^ have been considered real physical problems for which various 
technical approaches have been proposed. These problems are, for example, 
the galaxy-cluster mismatch, luminosity segregation, the richness-clustering 
relation and the linear non-linear evolution of the perturbations corresponding 
to the "small" or "large" amplitudes of fluctuations. We can now see that 
all this problematic situation is not real and it arises only from a statistical 
analysis based on inappropriate and too restrictive assumptions that do not 
find any correspondence in the physical reality. It is also important to note 
that, even if the galaxy distribution would eventually became homogeneous at 
larger scales, the use of the above statistical concepts is anyhow inappropriate 
for the range of scales in which the system shows fractal correlations as those 
shown in Fig.0. 



3 Radial density 

In the previous sections we have discussed the methods that allow one to mea- 
sure the conditional (average) density in real galaxy surveys. This statistical 
quantity is an average one, since it is determined by performing an average 
over all the points of the sample. We have discussed the robustness and the 
limits of such a measurement. We have pointed out that the estimate of the 
conditional density can be done up to a distance Ref / which is of the order of 
the radius of the maximum sphere fully contained in the sample volume. This 
is because the conditional density must be computed only in spherical shells. 
This condition puts a great limitation to the volume studied, especially in the 
case of deep and narrow surveys, for which the maximum depth Rg can be one 
order of magnitude, or more, than the effective depth Reff- 

We discuss here the measurement of the radial density in VL samples @. 
The determination of such a quantity allow us to extend the analysis of the 
space density well beyond the depth Reff- The price to pay is that such a 
measurement is strongly affected by finite size spurious fluctuations, becauseJi 
is not an average quantity. These finite size effects require a great cautionllZI: 
the behaviour of statistical quantities (like the radial density and the counts of 
galaxies as a function of the apparent magnitude) that are not averaged out, 
present new and subtle problems. 
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3. 1 Finite size effects and the behavior of the radial density 

In this section we discuss the general problem of the minimal sample size which 
is able to provide us with a statistically meaningful information. For example, 
the mass-length relation for a fractal, which defines the fractal dimension, is 

log(jV(< r)) 

D = hm — (28 

r^oo log(r) 

However this relation is properly defined only in the asymptotic limit, because 
only in this limit the fluctuations of fractal structures are self-averaging. A 
fractal distribution is characterized by large fluctuations at all scales and these 
fluctuations determine the statistical properties of the structure. If the struc- 
ture has a lower cut-off, as it is the case for any real fractal, one needs a 
"very large sample" in order to recover the statistical properties of the distri- 
bution itself. Indeed, in any real physical problem we would like to recover the 
asymptotic properties from the knowledge of a finite portion of a fractal and 
the problem is that a single finite realization of a random fractal is affected by 
finite size fluctuations due to the lower cut-off. 

In a homogeneous distribution we can define, in average, a chajjacteristics 
volume associated to each particle. This is the Voronoi volume Ii3 whose 
radius £v is of the order of the mean particle separation. It is clear that the 
statistical properties of the system can be defined only in volumes much larger 
than Vy. Up to this volume in fact we observe essentially nothing. Then one 
begins to include a few (strongly fluctuating) points, and finally, the correct 
scaling behavior is recovered (Fig.^). For a Poisson sample consisting of N 
particles inside a volume V then the Voronoi volume is of the order 

Vv-^^ (29) 

1/3 

and £y Vv .In the case of homogeneous distributions, where the fluctuations 
have small amplitude with respect to the average density, one readily recovers 
the statistical properties of the system at small distances, say, r > 5£v 

The case of fractal distribution is more subtle. Instead of the Voronoi 
length we can consider the average distance between nearest neighbors, but for 
clarity we proceed in the following way. Eqjl^ gives the mass-length relation 
for a fractal. In this case, the prefactor B is defined for spherical samples. If 
we have a sample consisting in a portion of a sphere characterized by a solid 
angle fl, we write 

N{< R) = Si?^^ . (30) 
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Figure 6: Behavior of the density computed from one point, in the case of a fractal distribu- 
tion. At small distances below the average mean separation between neighbor galaxies, one 
finds no galaxies. Then the number of galaxies starts to grow, but this regime is strongly 
affected by finite size fluctuations. Finally the correct scaling region r !v X is reached. In 
the intermediate region the density can be approximated roughly by a constant value. This 
leads to an apparent exponent D « 3. This exponent is not real but just due to the size 

effects. 
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In the case of a finite fractal structure, we have to take into account the 
statistical fluctuations. We can identify two basic kinds of fluctuations: the 
first ones are intrinsic f{R) and are due to the highly fluctuating nature of 
fractal distributions while the second ones, P{R), are Poissonian fluctuations. 
Concerning the first ones, one has to consider that the mass-length relation 
is a convolution of fluctuations which are present at all scales. For example 
one encounters, at any scale, a large scale structure and then a huge void: 
these fluctuations affect the power law behavior of N{< R). We can quantify 
these effects as a modulating term around the expected average given by Eq.pO|. 
Therefore, in the observations from a single point "i" we have 

[N{< i?)], = Bi?^^ • fn{R, Sn) . (31) 

In general it is more useful to focus on the behaviour of a local quantity (as 
the number of points in shells) rather than an integrated one. However for the 
purpose of the present discussion the approximation given by Eq.^ is rather 
good. This equation shows that the amplitude of A^(< r) is related to the 
amplitude of the intrinsic fluctuations and not only to the lower cut-off B. 
The amplitude of the modulating term is small, compared with the expected 
value of N{< R) 

V\IiW<BR''^. (32) 

In general this fluctuating term depends on the direction of observation Q and 
on the solid angle of the survey SH, so that f{R) = fQ{R,Sft). If we have a 
spherical sample we get 

^^/a(i?,5r!) = /4.(i?). (33) 



n 



In general we expect that fn{R, 6il) > fi.^{R), so that larger is the solid angle 
and smaller is the effect this term. If we perform the ensemble average of this 
fluctuating term we can smooth out its effect and we have then 

{k.{R))^ = 1 (34) 

where the average refers to all the occupied points in the sample. In 

such a way the conditional density, averaged over all the points of the sam- 
ple, has a single power law behavior. We stress that according to Eq.|l| the 
fluctuations in the integrated number of points in a fractal, are proportional 
to the number of points itself, raljhpr^jl:^^ to the root mean square as in a 



poissonian distribution. In generalli3E3tiil^E2l it is possible to characterize these 
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intrinsic fluctuations as log-periodic oscillations in the power law behavior. By 
performing an ensemble average as in Eq.^ these oscillations can be smoothed 
out. However for the purpose of the present paper, we limit our discussion to 
the approximation of Eq.^, without entering in more details. 

The second P{R) term is an additive one, and it takes into account spurious 
finite size fluctuations is simply due to shot noise. In this case we have that 

V\PiW >BR"^ if R<X (35) 
47r 

while P{R) ~ for R> X. The ensemble average is, again, expected to be 

{P{R}h = . (36) 

This term becomes negligible if the shot noise fluctuations are small: for ex- 
ample, if 

[iV(< i?)],; > 10^[N{< R)], . (37) 

From this condition and Eqjs^ we can have a condition on A (neglecting the 
effect of /(i?)): 

W-(^O^^Y (38) 



Bn J 

The minimal statistical length X is an explicit function of the prefactor B and of 
the solid angle of the survey fi. This scale is a lower limit for the scaling region 
of the distribution: the effect of intrinsic fluctuations, described by f{R), which 
are in general non negligible, can modulate the distance at which the scaling 
region is reached. This length depends also, but weakly, on the particular 
morphological features of the sample. Therefore it is important to stress that 
Eq.^ gives an order of magnitude for (A), where we intend the average value 
over all the possible directions of observations. In different directions one can 
have different values for (A), because of the effect of /(r). 

In the case of real galaxy catalogs we have to consider the luminosity 
selection effects. In a VL sample, characterized by an absolute magnitude 
limit Miim, the mass- length relation Eq.^, can be generalized as 

N{R, Mu^) = BR°^yj{Mu„,) (39) 

where tp{Miim) is the probability that a galaxy has an absolute magnitude 
brighter than Mum 

f*^'™ (t>{M)dM 
< ^{Mum) = Z . < 1 (40) 
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Survey 


n{sr) 


\{h-H.Ipc) 


CfAl 


1.8 


15 


CfA2 (North) 


1.3 


20 


SSRSl 


1.75 


15 


SSRS2 


1.13 


20 


PP 


1 


40 


LEDA 


2 TT 


10 


IRAS 


2 TT 


15 


ESP 


0.006 


300 



Table 1: The minimal statistical length A for several redshift surveys 



where (t>{M) is the Schecther luminosity function (normaHzed to unity) and 
^'(oo) is the normahzing factor 

*(oo) = / (t>{M)dM (41) 



where Mmin is the fainter absolute magnitude surveyed in the catalog (usually 

Mrmn « "10 --11). 

It is possible to compute the intrinsic prefactor B from the knowledge 
of the conditional density r(r) (Eq.^tla) computed in the VL samples and 
normalized for the luminosity factor (Eq.Eo). In the various VL subsamples of 
Perseus-Pisces, CfAl, and other redshift surveys we find that 

B « 10 ^ 15ih-Hlpc)-" (42) 

depending on the parameters of the Schecther function M* and S. From Eq.^, 
Eq.^ and Eq.^ we obtain for a typical volume limited sample with Mum ~ 

M*, 

(A) « (^»-^^")r'^'"^ . (43) 

This is the value of the minimal statistical length that we use in what follows. 
In Tab. ^ we report the value of A for several redshift surveys. While in 
the case of CfAl, SSRSl, PP, LEDA and ESP we have checked that there 
is a reasonable agreement with this prediction, the CfA2 and SSRS2 redshift 
surveys are not sill published and hence in these cases we can predict the value 
of A. 
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Given the previous discussion, we can now describe in a very simple way 
the behavior of [N{< r)]i, i.e. the mass length relation measured from a generic 
point "i". Given a sample with solid angle Q, we can approximate the effect 
of the intrinsic and shot noise fluctuations in the following way: 

(A^(< r)), = if r < A (44) 

i.e. the density is constant up to A, while 

(iV(< r)),; = Br^ if r > A (45) 
so that by the condition of continuity at A we have 

^1 - ^ (46) 

This simple approximation is very useful in the following discussion, especially 
for the number counts. 

To clarify the effects of the spatial inhomogeneities and finite effects we 
have studied the behavior of the galaxy radial (number) density in the VL 
samples, i.e. the behavior of (using Eq.^ and Eq.^ 

nvLir) = - ^Br^--'i^{Mum) ■ fn{r, Sn) (47) 

V(r) 47r 

One expects that, if the distribution is homogeneous the density is constant, 
while if it is fractal it decays as power law. 

When one computes the conditional average density, one indeed performs 
an average over all the points of the survey. In particular, as we have already 
discussed, we limit our analysis to a size defined by the radius of the maximum 
sphere fully contained in the sample volume, and we do not make use any 
treatment of the sample boundaries. On the contrary Eq.^ is computed only 
from a single point, the origin. This allows us to extend the study of the spatial 
distribution up to very deep scales. The price to pay is that this method is 
strongly affected by statistical fluctuations and finite size effects. 

The effect of the finite size spurious fluctuations for a fractal distribution 
is shown Fig.^: at small distances one finds almost no galaxies because we are 
below the average separation between neighbor galaxies £. Then the number 
of galaxies starts to grow, but this regime is strongly affected by finite size 
fluctuations. Finally the correct scaling region r « A is reached. This means, 
for example, that if one has a fractal distribution, there is flrst a raise of 
the density, due to finite size effects and characterized by strong fluctuations, 
because no galaxies are present before a certain characteristic scale. Once one 
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enters in the correct scaling regime for a fractal the density starts to decay as 
a power law. So in this regime of raise and fall with strong fluctuations there 
is a region in which the density can be approximated roughly by a constant 
value. This leads to an apparent exponent 13 « 3, so that the integrated 
number grows as N{< r) ~ r'^. This exponent is therefore not a real one 
but just due to finite size fluctuations. Of course, depending on the survey 
orientation in the sky, one can get an exponent larger or smaller than 3, but 
in general this is the more frequent situation (see below). Only when a well 
defined statistical scaling regime has been reached, i.e. for r > A, one can 
find the genuine scaling properties of the structure, otherwise the behavior is 
completely dominated by spurious finite size effects (for seek of clarity in this 
discussion we do not consider the effect of f(r) in Eq.^l|) 

The question of the difference between the integration from the origin (ra- 
dial density) and correlation properties averaged over all points lead us to a 
subtle problem of asymmetric fluctuations in a fractal structure. From our 
discussion, exemplified by FigJ|, the region between i and A corresponds to 
an underdensity with respect to the real one. However we have also showed 
that for the full correlation averaged over all the points, as measured by r(r), 
the correct scaling is recovered at distances appreciably smaller than A. This 
means that in some points one should observe an overdensity between i and 
A. However, given the intrinsic asymmetry between filled and empty regions 
in a fractal, only very few points show the overdensity (a fractal structure is 
asymptotically dominated by voids). These few points nevertheless have, in- 
deed, an important effect on the average values of the correlations. This means 
that, in practice, a typical points shows an underdensity up to A as shown in 
Fig.g. The full averages instead converge at much shorter distances. This 
discussion shows the peculiar and asymmetric nature of finite size fluctuations 
in fractals compared to the symmetric Poissonian case. For homogeneous 
distributioiO the situation is in fact quite different. Below the Voronoi length 
iy there are finite size fluctuations, but for distances r > (2 4)£^ the cor- 
rect scaling regime is readily found for the density. In this case the finite size 
effects do not affect too much the properties of the system because a Poisson 
distribution is characterized by small amplitude fluctuations. 

As an example, we show the behaviour of the radial density in the Perseus- 
Pisces redshift survey (PPRS). We have computed the n(r) in the various 
VL samples of PPRS, and we show the results in Fig.Q. In the less deeper 
VL samples (VL70, VL90) the density does not show any smooth behavior 
because in this case the finite size effects dominate the behavior as the distances 
involved are r < A (Eq.^3|) . At about the same scales we find a very well defined 
power law behavior by the correlation function analysis. In the deeper VL 
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PERSEUS-PISCES 




Figure 7: The spatial density n(r) computed in the VL sample cut at 70, 90, 110, 130h~^ Mpc 
. In the case of VL70 and VL90 the density is dominated by large fluctuations and it has 
not reached the scaling regime. In the samples VLllO and VL130 the density is dominated 
by large fluctuations only at small distances, while at larger distances, after the Perseus 
Pisces chain at 50h~^ Mpc, a very well defined power law behavior is shown, with the same 
exponent of the correlation function (i.e. D = 2) 



samples (VLllO, VL130) a smooth behavior is reached for distances larger than 
the scaling distance {fl — 0.9 sr) r « A ~ bOh'^^Mpc. The fractal dimension 
is Z? « 2 as one measures by the correlation function. For relatively small 
volumes it is possible to recover the correct scaling behavior for scales of order 
of £ (instead of ~ 10^) by averaging over several samples or, as it happens in 
real cases, over several points of the same sample when this is possible. Indeed, 
when we compute the correlation function we perform an average over all the 
points of the system even if the VL sample is not deep enough to satisfy the 
condition expressed by JEq.^. In this case the lower cut-off introduces a limit 
in the sample statistics □. 

The case of LCRS and ESP are discussed in detail in Sylos Labini et 
al.(1997) and the results are shown in Figj^ 

3.2 Pencil beams 

Deep "pencil beams" surveys cover in general very narrow angular size ('^ 
1°) and extend to very deep depths (z > 0.2). These narrow shots through 
deep space provide a confirmation of strong inhomogeneities in the galaxy 
distribution. 

One of the most discussed results obtained from pencil-beams surveys has 
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been the claimed detection of a typical scale in the distribution of gaJaxy struc- 
tures, corresponding to a characteristic separation of 128h~'^Mpcc3. However 
in the last five years several other surveys, in different regions of the sky, do 
not find auM evidence for such a periodicity. In particular, Bellanger & De 
LapparentcH, by analyzing a sample of 353 galaxies in the redshift interval 
0.1 < z < 0.5, concluded that these new data contain no evidence for a pe- 
riodic signal on-a scale of 128h~^Mpc. Moreover they argue that the low 
sampling rate ofo is insufficient for mapping the detailed large-scale structure 
and it is the real origin of the apparent periodicity. 

On the other hand, Willmer et-aZ.(1994) detect four of the five nearest 
peaks oLthe galaxies detected byE3, because their survey is contiguous to 
that ofc2l, in the sky region near the north Galactic pole. Moreover Ettori et 
aL(1996) in a survey oriented in three small regions around the South Galactic 
Pole do not find ami statistically periodic signal distinguishable from noise. 
Finally Cohen et altB by analyzing a sample of 140 objects up to z 0.8, find 
that there is no evidence for periodicity in the peak redshifts. 

In a pencil-beam survey one can study the behavior of the linear den- 
sity along a tiny but very long cylinder. The observed galaxy distribution 
corresponds therefore to the intersection of ±he full three dimensional galaxy 
distribution with one dimpnsional cylinder EJ. In this case, from the law of 
codimension additivity Eilu, one obtains that the fractal dimension of the in- 
tersection is given by 

Di = D + dpc-dKD + l- 3aiO (48) 

where D is the galaxy distribution fractal dimension, embedded in a d = 3 
Euclidean space, and dpc = 1 is the dimension of the pencil beam survey. This 
means that the set of points visible in a randomly oriented cylinder has dimen- 
sion Di « 0. In such a situation the power law behavior is no longer present 
and the data should show a chaotic, featureless nature strongly dependent on 
the beam orientation. If the galaxy distribution becomes instead homogeneous 
above some length, shorter than the pencil beam depth, one has the regular 
situation Di — 1 and a well defined density must be observed. 

We stress that in any of the available pencil beams surveys, one can de- 
tect tendency towards an homogeneous distribution. Rather, all these surveys 
show a very fluctuating signal, characterized by the presence of galaxy struc- 
tures. Some authorsllj claim to detect the end greatness (i.e. that the galaxy 
structures in the deep pencil beams are not, so different from those seen in 
nearby sample - as the Great Wall), or thattZl the dimension of voids does not 
scale with sample size, by the visual inspection of these surveys. However one 
should consider in these morphological analyses, that one is just looking at a 
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convolution of the survey geometry, which in general are characterized by very 

narrow solid angles, and large scale structures, and that in such a situation, a 
part from very favorable cases, one may detect portions of galaxy structures 
(or voids). 

Finally wc note that if the periodicity would be present, the amplitude of 
the different peaks is very different from each other, and an eventual transition 
to homogeneity in a periodic lattice, should be, for example, ten times the 
lattice parameter, i.e. > 1000h~^Mpc ! 

4 Number Counts 

The most complete information about galaxy distribution comes from the full 
three dimensional samples, while the angular catalogs have a poorer qualita- 
tive information, even if usually they contain a much larger number of galaxies. 
However, one of the most important tasks in observational astrophysics, is the 
determination of the log — log S relation for different kind of objects: galax- 
ies in the various spectral band (from ultraviolet to infrared), radio-galaxies. 
Quasars, X-ray sources and 7-ray bursts. This relation gives the number N{S) 
(integral or differential) of objects, for unit solid angle, with apparent flux 
(larger than a certain limit) S. The determination of such a quantity avoids 
the measurements of the distance, which is always a very complex task. How- 
ever wc show in the following that the behavior of the log N — log S is strongly 
biased by some statistical finite size effects due to small scale fluctuations. 

The counts of galaxies as a function of the apparent magnitude are deter- 
mined from the Earth only, and hence it is not possible to make an average over 
different observers. As we have already discussed in the previous section, this 
kind of measurement is affected by intrinsic fluctuations that are not smoothed 
out at any scale. Moreover at small scale there are finite size effects which may 
seriously perturb the behavior of the observed counts. Following the simple 
argument about the behaviour of the radial density we have presented in the 
previous section, we consider here the problem of the galaxy-number counts. 

We present in this section a new interpretation of this basic relation at the 
light of the highly inhomogeneity nature of galaxy distribution, and we show 
its compatibility with the behavior of counts of galaxies in different frequency 
bands, radiogalaxies. Quasars and X-ray sources. Our conclusion is that the 
counts of all these different kind of objects are compatible with a fractal distri- 
bution of visible matter up to the deepest observable scale. 
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Galaxy number counts data 

Historically SEl the oldest type of data about galaxy distribution is given 
by the relation between the number of observed galaxies A^(> S) and their 
apparent brightness S. It is easy to show that, under very general assumptions 



one gets (see Sec. 4.2) 



A^(> S) - S-- (49) 

where D is the fractal dimension of the galaxy distribution. Usually this re- 
lation is written in terms of the apparent magnitude m {S ^ 2^Q-o.4?n _ j^Q^g 
that bright galaxies correspond to small to). In terms of to, Eq.^ becomes 

log iV(< to) ~ Q!TO (50) 

with a = Z3/5&Ei. Note that a is the coefficient of the exponential behavior 
of Eq.^ and we call it " exponent" even though it should not be confused with 
the exponents of power law behaviors. In Fig.|^ we have collected jjlliith'j' IpYfflt 
observations of N{< to) versus to in the i3-spectral-band TOBE3'E3Eill23SE3t3. 
At bright and intermediate magnitudes (12 < niB ^ 18), corresponding to 
small redshift { z < 0.2), one obtains a w 0.6, while from tub ^ 19 up 
to iriB ^ 28 the counts axe, welt-fitted by a smaller exponent with a w 0.4. 
The usual interpretation liS'HEjfiJ'a is that a ~ 0.6 corresponds to « 3 
consistent with homogeneity, while at large scales galaxy evolution and space 
time expansion effects are invoked to explain the lower value a « 0.4. On 
the basis of the previous discussion of the VL samples, this interpretation is 
untenable. In fact, we know for sure that, at least up to ~ 150h~^Mpc there 
are fractal correlations, as we have discussed in the previous sections, so one 
would eventually expect the opposite behavior. Namely small value of a ~ 0.4 
(consistent with Z? w 2) at small scales followed by a crossover to an eventual 
homogeneous distribution at large scales {a « 0.6 and Z? « 3). 

The situation is therefore similar in the different spectral bands. The puz- 
zling behavior of the GNC represents an important apparent contradiction we 
find in the data analysis. We argue here that this apparently contradictory 
experimental situation can be fully understood on the light of the small scale 
effects in the space distribution of galaxies. For example a fractal distribution 
is non analytic in every occupied point: it is not possible to define a meaningful 
average density because we are dealing with intrinsic fluctuations which grow 
with as the scale of the system itself. This situation is qualitatively differ- 
ent from an homogeneous picture, in which a well defined density exists, and 
the fluctuations represent only small amplitude perturbations. The nature of 
the fluctuations in these two cases is completely different, and for fractals the 
fluctuations themselves define all the statistical properties of the distribution. 
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Figure 8: The galaxy number counts in ttie B-band, from several surveys. In the range 
12 < m < 19 the counts show an exponent a 0.6±0.1, while in the range 19 ^ ^ 28 the 
exponent is o k:^ 0.4. The solid line is computed from the determination of the amplitude of 
the conditional density at small scale, the fractal dimension D = 2.2, and from the knowledge 
of the luminosity function. The distance is computed for a galaxy with M = —16 and we 
have used Hq = 75kmsec~^ Mpc~^ . 
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This concept has dramatic consequences in the fohowing discussion as weh as 
in the determination of various observable quantities, such as the amphtude of 
the two point angular correlation function. It is worth to notice that the small 
scale effects are usually neglected in the study of fractal structures because one 
can generate large enough (artificial) structures to avoid these problems. In 
Astrophysics the data are instead intrinsically limited and a detailed analysis 
of finite size effects is very important. We discuss the problems of finite size ef- 
fects in the determination of the asymptotic properties of fractal distributions, 
considering explicitly the problems induced by the lower cut-off. 

4-2 Galaxy counts: basic relations 

We briefly introduce some basic relations which are useful later. We can start 
by computing the expected GNC in the simplest case of a magnitude limited 
(ML) sample. A ML sample is obtained by measuring all the galaxies with 
apparent magnitude brighter than a certain limit mum. In this case we have 
(for m < mii,n) 

N{< m) =^ S$((X))10T™ (51) 

where 

/oo 
(?i.(Af)10-^(*^+25)dM . (52) 
-oo 

We consider now the case of a volume limited (VL) sample. A VL sample 
consists of every galaxy in the volume which is more luminous than a certain 
limit, so that in such a sample there is no incompleteness for an observational 
luminosity selection effectQ'EJ. Such a sample is defined by a certain maximum 
distance R and an absolute magnitude limit given by: 

Mii„i = TOHm - 5 logio i? - 25 - A{z) (53) 

(nT-iim is the survey apparent magnitude limit). By performing the calculations 
for the number-magnitude relation, we obtain 

iV(< m) = A{m) ■ 10"?™ + C(m) (54) 

where v4(m) is 

A{m)^B 0(M)lO--(A^+25)^^ (55) 

JMim) 

and M{m) is given by M(rn) — m — 51og(_R) — 25, and it is a function of m. 
The second term is 

/M(m) 
(j){M)dM . (56) 
-OO 
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This term, as A{m), depends on the VL sample considered. We assume a 
luminosity function with a Schecther shape. For M{m) > M* we have that 
C(m) w 0, and A(m) is nearly constant with m. This happens in particular 
for the deeper VL samples for which Mum ~ M* . For the less deeper VL 
{Miim, > M*) samples these terms can be considered as a deviation from a 
power law behavior only for m mum- If one has 4>{M) = 5[M — Mq) then 
it is simple to show that log(A^(< m)) ^ {D/5)m also in each VL sample. 

4-3 Galaxy counts in redshift surveys 

We have studied the GNC in the Perseus-Pisces redshift survey @ in order to 
clarify the role of spatial inhomogeneities and finite size effects. In the previous 
sections we have analyzed the spatial properties of galaxy distribution in this 
sample by measuring the conditional (average) density and the radial density. 
Let us briefly summarize our main results. 

When one computes the conditional average density, one indeed performs 
an average over all the points of the survey. On the contrary the radial density 
is computed only from a single point, the origin. This allows us to extend the 
study of the spatial distribution up to very deep scales: the price to pay is 
that this method is strongly affected by statistical fluctuations and finite size 
effects. Analogously, when one computes A'^(< m), one does not perform an 
average but just counts the points from the origin. As in the case of the radial 
density n{r) also A^(< to) is strongly affected by statistical fluctuations due to 
finite size effects, as well as intrinsic oscillations that are not smoothed out. 
We are now able to clarify how the behavior of A^(< m), and in particular its 
exponent, are influenced by these effects. 

We show in Fig.^ and Fig.^ the behavior of N(< to) respectively for the 
various VL samples and for the whole magnitude limit sample. In VL60 there 
are very strong inhomogeneities in the behavior of n{r) (see Fig.|^) and these 
are associated with a slope a « 0.6 for the GNC. (The flattening for to mum 
is just due to a luminosity selection effect that is explained in Sec. [4.2[) . For 
VLIIO the behavior of the density is much more regular and smooth, so that 
it shows indeed a clear power law behavior. Correspondingly the behavior of 
iV(< to) is well fitted by a w 0.4. Finally the whole magnitude limit sample 
is again described by an exponent a « 0.6. 

We have now enough elements to describe the behavior of the GNC. The 
first point is that the exponent of the GNC is strongly related to the space 
distribution. Indeed- what has never been taken into account before is the role 
of finite size effectstZl. The behavior of the GNC is determined by a convolution 
of the space density and the luminosity function, and the space density enters 
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Figure 9: (a) The Number counts A''(< m) for the VL sample VL60. The slope is a 0.6. 
In this case occurs a flattening for m »Ti;im- (b) The Number counts A''(< m) for the 
VL sample VLllO. The slope is a ^ 0.4, a part from the initial fast growth due to weak 
statistics. This behavior corresponds to a well defined define power law behavior of the 
density with exponent D ss 5a » 2. 
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Figure 10: The Number counts A'^(< m) for the whole magnitude limit sample. The slope 
is a Si 0.6 and it is clearly associated only to fluctuations in the spatial distribution rather 
than to a real homogeneity in space. 

in the GNC as an integrated quantity. The problem is to consider the correct 
space density in the interpretation of data analysis. In fact, if the density has a 
very fluctuating behavior in a certain region of length scale, as in the case shown 
in Fig.0, its integral over this range of length scales is almost equivalent to a 
flat one. This can be seen also in Fig|6| at small distances one finds almost no 
galaxies because one is below the mean minimum separation between neighbor 
galaxies. Then the number of galaxies starts to grow, but this regime is strongly 
afli'ected by finite size fluctuations. Finally the correct scaling region r w A is 
reached. This means for example that if one has a fractal distribution, there is 
flrst a raise of the density, due to flnite size effects and characterized by strong 
fluctuations, because no galaxies are present before a certain characteristic 
scale. Once one enters in the correct scaling regime for a fractal the density 
becomes to decay as a power law. So in this regime of raise and fall with strong 
fluctuations there is a region in which the density can be approximated roughly 
by a constant value. This leads to an apparent exponent Z? w 3, so that the 
integrated number grows as A^(< r) ~ and it is associated in terms of GNC, 
to a « 0.6 (Fig.pr]). This exponent is therefore not a real one but just due 
to finite size fiuctuations. Only when a well defined statistical scaling regime 
has been reached, i.e. for r > A, one can find the genuine scaling properties 
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Figure 11: The number counts N(< m) together with the behavior of the space density. At 
small scale the density is characterized by having strong fluctuations which lead to a slope 
a f» 0.6. This is clearly associated only to fluctuations in the spatial distribution rather than 
to a real homogeneity in space. At larger scales (faint end) the correct scaling behavior is 

recovered and a = D/5. 
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of the structure, otherwise the behavior is completely dominated by spurious 
finite size effects. In the VL samples where n{r) scales with the asymptotic 
properties (Fig.^ the GNC grows also with the right exponent {a = D/5). 

If we now consider instead the behavior of the GNC in the whole magnitude 
limit (hereafter ML) survey, we find that the exponent is a « 0.6 (Fig.p^). This 
behavior can be understood by considering that at small distances, well inside 
the distance A defined by Eq.|3[ the number of galaxies present in the sample 
is large because there are galaxies of all magnitudes. Hence the majority of 
galaxies correspond to small distances (r < A) and the distribution has not 
reached the scaling regime in which the statistical self-averaging properties 
of the system are present. For this reason in the ML sample the finite size 
fluctuations dominate completely the behavior of the GNC. Therefore this 
behavior in the ML sample is associated with spurious finite size effects rather 
than to real homogeneity. We discuss in a more quantitative way the behavior 
in ML surveys later. 

4-.4 Test on finite size effects: the average N{< m) 

To prove that the behavior found in Fig.p^, i.e. that the exponent a ~ 0.6 is 
connected to large fluctuations due to finite size effects in the space distribution 
and not to real homogeneity, we have done the following test. We have adopted 
the same procedure used for the computation of the correlation function, i.e. 
we make an average for iV(< m) from all the points of the sample rather than 
counting it from the origin only. 

To this aim we have considered a VL sample with N galaxies and we have 
built N independent flux-limited surveys in the following way. We consider 
each galaxy in the sample as the observer, and for each observer we have 
computed the apparent magnitudes of all the other galaxies. To avoid any 
selection effect we consider only the galaxies which lie inside a well defined 
volume around the observer. This volume is defined by the maximum sphere 
fully contained in the sample volume with the observer as a center. 

Moreover we have another selection effect due to the fact that our VL sam- 
ple has been built from a ML survey done with respect to the origin. To avoid 
this incompleteness we have assigned to each galaxy a constant magnitude M . 
In fact, our aim is to show that the inhomogeneity in the space distribution 
plays the fundamental role that determines the shape of the iV(< m) relation, 
and the functional form of the luminosity function enters in Eq.|5^ only as an 
overall normalizing factor. 

Once we have computed Ni{< m) from all the points z = 1, .., iV we then 
compute the average. We show in Fig.O the results for VL60 and VLllO: a 
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Figure 12: Left panel The average N(< m) in the VL sample VL60. The squares crosses refer 
to the average N{< m) computed assigning to all the galaxies the same absolute magnitude 
Mo = M*. The reference line has a slope a = 0.4. Right panel The average N(< m) in the 
VL sample VLllO. The squares crosses refer to the average Af(< m) computed assigning 
to all the galaxies the same absolute magnitude Mq = M* . The reference line has a slope 

a = 0.4 

very well defined exponent a = D/b ~ 0.4 is found in both cases. This is in 
fully agreement with the average space density (the conditional average density 
r(r)) that shows D « 2 in these VL samples. 

4-5 Behavior of galaxy counts in magnitude limited samples 

We are now able to clarify the problem of ML catalogs. Suppose to have 
a certain survey characterized by a solid angle and we ask the following 
question: up to which apparent magnitude limit mum we have to push our 
observations to obtain that the majority of the galaxies lie in the statistically 
significant region (r > A) defined by Eq.^. Beyond this value of mum we 
should recover the genuine properties of the sample because, as we have enough 
statistics, the finite size effects self-average. From the previous condition for 
each solid angle f2 we can find an apparent magnitude limit m/im so that finally 
we are able to obtain = i^Umi^) in the following way. 

In order to clarify the situation, we can now compute the expected value 
of the counts if we use the approximation for the behavior of the mass-length 
relation given by Eqs. ^^ji^ . Suppose, for seek of clarity, also that (j>{M) = 
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5{M - Mo), with Mo = -19. We define 

A = 100-2(™a-J>^-25) (57) 

where A is given by Eq.O. Then the differential counts are given by 



dm J ^ 5 



TOA (58) 



and 



dN\ ^ log^^^ . ^0-™ . if ^ > (59) 

dm I ■ 5 



If Mo = -19 and A ^ j^ih^^Mpc) we have 

5 

wiim = m(ri) « 14 — — logil (60) 

In order to give an estimate of such an effect if <I>(M) has a Schechter like 
shape, we can impose the condition that, in a ML sample, the peak of the 
selection function, which occurs at distance rpeak ^ satisfies the condition 

Tpeak > A (61) 

where A in the minimal statistical length defined by Eq.^. The peak of the 
survey selection function occurs for M* « —19 and then we have rpeak ~ 
^Q_LLm ^ From the previous relation and from Eq|6^ and Eq.^ we easily 
recover Eqg 

The magnitude m(fl) separates the 0.6 behavior, strongly dominated by 
intrinsic and shot noise fiuctuations, from the asymptotic 0.4 behavior. Of 
course this is a crude approximation in view of the fact that A has not a well 
defined value, but it depends on the direction of observation and not only on 
the solid angle of the survey. However the previous equations gives a reasonable 
description of real data. 

We show in Fig.^ the condition given by Eq.^. From the previous figure 
it follows that for m > 19 the statistically significant region is reached for 
almost any reasonable value of the survey solid angle. This implies that in deep 
surveys, if we have enough statistics, we readily find the right behavior {a = 
D/5) while it does not happens in a self-averaging way for the nearby samples. 
Hence the exponent a ~ 0.4 found in the deep surveys (m > 19) is a genuine 
feature of galaxy distribution^ and corresponds to real correlation properties. 
In the nearby surveys m < 17 we do not find the scaling region in the ML 
sample for almost any reasonable value of the solid angle. Correspondingly 
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Figure 13: If a survey defined by the ajjparent magnitude limit rnii„^ and the solid angle Q 
lie in the statistically significant region it is possible to obtain the self-averaging properties 
of the distribution also witii tiie integral from the vertex. Otherwise one needs a redshift 
survey which contains the three dimensional information, and then one can perform average. 
Only in this way it is possible to smooth out the finite size effects. 



35 



the value of the exponent is subject to the finite size effects, and to recover the 
real statistical properties of the distribution one has to perform an average. 

From the previous discussion it appears now clear why a change of slope 
is found at m ~ 19: this is just a reflection of the lower cut-off of the fractal 
structure and in the surveys with m;,;,„ > 19 the self-averaging properties of the 
distribution cancel out the finite size effects. This result depend very weakly 
on the fractal dimension D and on the parameters of the luminosity function 
6 and M* used. Our conclusion is therefore that the exponent a « 0.4 for 
m > 19 is a genuine feature of the galaxy distribution and it is related to a 
fractal dimension D k, 2, which is found for m < 19 in redshift surveys only 
by performing averages. We note that this result is based on the assumption 
that the Schecther luminosity function holds also at high redshift, or, aL least 
to m ~ 20. This result is confirmed by the analysis of Vettolani et a£j who 
found that the luminosity fymction up to z ^ 0.2 is in excellent agreement with 
that found in local surveysu. 

Finally a comment on the amplitude of counts. In Fig.|^ the solid line 
represents the behavior of Eqjs^. The prefactor B has been determined by the 
full correlation analysis while the fractal dimension is D = 2.2. The parameter 
of the luminosity function are 5 = —1.1 and M* = —19.5 as usual. The 
agreement at faint magnitudes (m > 19) is quite good. At bright magnitudes 
one usually underestimates the number of galaxies because of finite size effects, 
even in some cases the number counts can be larger than the predicted value. 
This is related to the asymmetry of the fluctuations in a fractal structure. For 
example we report in Fig..^ various determinations of the GNC in different 
regions of the sky (fromtill) . The fact that at faint magnitudes the behavior is 
quite regular, is due to the smoothing of spatial fluctuations for the luminosity 
effect. Namely at a given apparent magnitude there are contribution from 
galaxies located at very different distance, as the luminosity function of galaxies 
is spread over several decades of luminosities. On the other hand at the bright 
end there are contributions only from nearby galaxies, and is such a way the 
space finite size fluctuations are not smoothed out. 

In the data shown in Fig.^ K-corrections have been not applied. Such 
corrections must be present because of the Hubble distance-redshift relation. 
Namely the spectrum of a certain galaxy at redshift z is shifted towards red 
of a certain amount, according to the Hubble law. There are several galaxies 
(E/SO) with steep spectra and hence for these one detects a lower value of the 
intensity of the "true" one because of the shift of the maximum of the spectrum. 
However there are several other galaxies (Sdm, Scd) with flat spectra. In some 
case the shift due to the Hubble law may produce an higher intensity while in 
other a decrease of the apparent flux. The K-corrections are model dependent 
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Figure 14: Determinations of the GNC in various sky regions (from Bertin & Dennefeld, 
1996). The bright end wildly fluctuates around the predicted value, while at fainter magni- 
tudes the behavior is quite smooth and is the same for the different fields. 



corrections: the sign can be in both directions, i.e. towards an increasing or a 
decreasing of the absolute magnitudecj. Moreover the effect of such corrections 
gin general not so important for the distortion of the number counts relation 

In observational astrophysics there are a lot of data which are only angular 
ones, as the measurements of distances is general a very complex task. We 
briefly discuss here the distribution of radiogalaxies, Quasars and the 7-ray 
burst distribution. Our conclusion is that all these data are compatible with 
a fractal structure with D w 1.6 2.0 similar to that of galaxies. 



4-. 6 Radio galaxies 

As an example of counts of objects, we consider the case of radio galaxies. The 
majority of catalog radio sources are extragalactic and some of the strongest 
are at cosmological distances. One of the most important information on radio 
galaxies distribution_has been obtained from the sources counts as a function 
of the apparent fluxCJ. Extensive surveys of sources have been made at various 
frequencies in the range 0.4 —IlGHz. In Fig.^ we show a compendium of 
sources counts at v — lAGhzcB. The differential counts are plotted against 
the apparent flux. In the bright flux region there is a deviation from a power 
law function, while for four decades the agreement with a fractal distribution 
with D « 1.8 is excellent. Such a behavior has been usually explained in 
literature as an effect of sources/space-time evolution. Here we propose that 
the radio galaxies are fractally distributed, as galaxies, with almost the same 
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Figure 15: Normalized differential sources counts at u = lAGhz. Abscissa log flux density 
(Jy). Ordinate log differential number of sources n{S). The solid line represents the behavior 
of a fractal structure with D = 1.8. The amplitude has been computed from the knowledge 
of the luminosity function (see text). The agreement is excellent, except in the bright fluxes 
region due to the presence of finite size fluctuations. (Prom Condon, 1984) . 
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fractal dimension. The deviation at bright fluxes in this picture is explained 
as a spurious effect due to the small scale fluctuations, as in the case of galcixy 
counts. In the other frequency bands the situation is nearly the same c3. 
The simple picture of a fractal distribution of radio sources is therefore fully 
compatible with the experimental situation. 

It is simple to show that the differential number of radio galaxies for unit 
flux (in Jansky) and unit steradian, in an Euclidean space, is given by: 

where we have taken, as usual, n{r) — {B^D / A-K)r^^^ and S = L/(47rr^) 
(L is the intrinsic luniinosity) . From the knowledge of the parameters of the 
luminosity functioned we have that a good approximation of the data reported 
in Fig.^ is given hy D = 1.8 and Bb. = 0.1Mpc~^ . Such a low value for 
Bii implies that the density of radio galaxies is about 100 lower than that of 
optical galaxies. When it will be possible to have a complete redshift sample 
of Radio galaxies, i.e. a well defined volume limited sample, one will be able 
to measure directly Br and D from the knowledge of the conditional density 
and hence to compare those values with the ones we have measured from the 
number flux relation. 

4-7 Compendium of counts 

In Tab.| we summarize the behavior of the number counts for various astro- 
physical objects. The small scale exponent corresponds to the bright end of 
the number counts relation, while the large scale exponent to the faint end. 
The " small scale exponent" is not a real exponent and wildly fluctuates for the 
different objects. On the contrary the large scale exponent is a real exponent 
and its value is in the range 1.8 < D < 2.2 for almost all the cases. The counts 
of all these different objects is therefore compatible with a fractal distribution 
in space. 

This implies a completely new interpretation of the counts. At small scale 
the exponent I? « 3 (i.e a « 0.6 in the case of magnitude counts) is due to 
finite size effects and not to a real homogeneity: this has be shown to be the 
case with very specific tests. On the other hand, at larger scale the value 
D K, 2 (i.e. a ~ 0.4) corresponds to the correct correlation properties of the 
samples, that we can find by the complete correlation analysis in the three 
dimensional space. This implies that galaxy evolution, modifications of the 
Euclidean Geometry and the K-corrections are not very relevant in the range 
of the present data. In addition the fact that the exponent a ~ 0.4 holds up to 



39 



Objects 


Small scale behavior 




Large scale behavior 




U-band' 


1<U<18 


D 




18 < [/ < 24 




i 2.5 


B-band^ 


12 < B < 18 


D 


5i 3 


18 < B < 28 


D 


« 2 


V-band^ 


^ 20 


D 


PS? 


22 < y < 25 


D « 


1.95 


R-band* 


15 < i? < 18 


D R 


i2.8 


20 < < 26 


D « 


1.85 


I-band^ 


? ^ / ^ 19 


D 


w? 


19 < / < 25 


D R 


i 1.7 


K-band'5 


12<i^ < 17 


D ^ 


3.35 


20 < if < 24 


D R 


i 1.6 


Quasars'' 


14.75 <B< 18.75 


D K 


i 4.4 


19 < B < 23 


D R 


i 1.5 


Radio^ V = lAGhz 


1 < ^ < 10 


D « 


3.44 


10-^ <s<\ 


D R 


i 1.8 


Radio^ u = OMGhz 


1 ^S'< 10 


D R 


i3.8 


10-3 < 5 < 1 


£> R 


i 1.6 


Radio^ = OAOSGhz 


1 ^ S* ^ 10 


D R: 


i 3.7 


10-3 <s<\ 


D R 


i 1.5 


Radio^ u = 5.0Ghz 


1 < ^ < 10 


D R 


J 3.4 


10-^ <s<\ 


D R 


i 1.8 


X-ray sources^ 


5 • 10-13 < 5 < 10-12 


D R 


i3.4 


10-16 < 5 < 5 • 10-13 


£> R 


i 1.8 


7 — ray bursts 


10 < 5 < 100 


D 


« 3 


10-1 < 5 < 10 


D R 


i 1.7 



Table 2: The exponents of counts for different kinds of astrophysical objects (see text). 
(Ref.1-6: see Table 8; Ref.7: Hartwick & Schade, 1990; Ref.8: Condon,1984; Ref.9: Hasinger 
et at, 1993; Ref.lO: Meegan et aZ.1995) 



magnitude 27 28 for galaxies seems to indicate that the fractal may continue 
up to distances ^ 4000/i-iMpc. Quite a remarkable fact if one considers that 
the Hubble radius of the Universe is supposed to be 4000/i-iMpc. Moreover 
such a behavior can be found for galaxies in the different photometric bands, 
as well as for other astrophysical objects. 

5 Conclusion 

The highly irregular galaxy distributions with large structures and voids strongly 
point to a new statistical approach in which the existence of a well defined av- 
erage density is not assumed a priori and the possibility of non analytical 
properties should be addressed specifically. The new approach for the study of 
galaxy correlations in all the available catalogues shows that their properties 
are actually compatible with each other and they are statistically valid samples. 
The severe discrepancies between different catalogues that have led various au- 
thors to consider these catalogues as not fair, were due to the inappropriate 
methods of analysis. 
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The correct two point correlation analysis shows well defined fractal cor- 
relations up to the present observational limits, from 1 to 1000h~^Mpc with 
fractal dimension D ~ 2. Of course the statistical quality and solidity of the 
results is stronger up to 100 2QQh~^Mpc and weaker for larger scales due 
to the limited data. It is remarkable, however, that at these larger scales one 
observes exactly the continuation of the correlation properties of the small and 
intermediate scales. 

From the theoretical point of view the fact that we have a situation char- 
acterized by self-similar structures implies that we should not use concept like 
^(r), ro, 5N/N and certain properties of the power spectrum, because they are 
not suitable to represent the real properties of the observed structures. In this 
respect also the N-body simulations should be considered from a new perspec- 
tive. One cannot talk about "small" or "large" amplitudes for a self-similar 
structure because of the lack of a reference value like the average density. The 
Physics should shift from "amplitudes" towards "exponent" and the methods 
of modern statistical Physics should be adopted. This requires the development 
of constructive interactions between two fields. 
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